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Abstract — The aim of this paper is to present a new fast- 
convergent numerically stable space-time adaptive processing 
(STAP) algorithm derived using a novel technique of feedback 
orthogonalization. The main advantages of this approach lie in 
its perfected stability to computational errors and faults which 
makes its real-time implementation on substantially faster and 
cheaper regular fixed-point processors possible. 

I. Introduction 

Most STAP algorithms designed for effective on-line 
application in fields such as radar, sonar, communications, 
seismology, radioastronomy, and medical imaging are 
required to have an extremely rapid rate of actual convergence. 
One way to speed up the real-time adaptive processing is to 
take advantage of faster cost-effective fixed-point parallel 
computing systems. 

Due to the fact that the optimum weights associated with 
such criteria as the LMS, maximum SINR, and MVDR differ 
from the optimum Wiener solution only by a scale gain factor 
and have equal output SINR's, we consider an adaptive space- 
time filter optimized to maximize the SINR at its output [1]. 

In terms of computational structure, the quadratic 
optimisation STAP algorithms can be broadly classified under 
asymptotically convergent iterative and finite convergent 
direct categories. Since the iterative LMS-type algorithms 
implementing local searching optimization methods are 
unacceptably slow in most practical scenarios [2], we focus 
our attention on the finite convergent algorithms based on 
direct methods of linear algebra. 

With respect to the number of iterations a numerical 
method requires to achieve an acceptable approximation to the 
optimum solution, the direct algorithms based on various 
types of matrix decomposition and recurrent inversion are 
considered to be very fast because under conditions of infinite 
precision arithmetic they converge in 2K iterations, where K is 
the filter dimension [3]. However, numerous results of 
convergence analyses [2], [4-7] have shown that the direct 
algorithms are explosively divergent in regular finite precision 
environments necessitating their implementation on slower 
expensive floating-point processors. 

The underlying cause of numerical instability of most direct 
STAP algorithms is in the strictly sequential feedback-less 
processing structure that results in uncontrollable propagation 
and accumulation of numerical errors in finite precision 
computing systems. 



The main goal of this paper is to present a new fast- 
convergent algorithm for real-time STAP designed to be 
perfectly robust in terms of processing errors and faults that 
provides its real-time implementation on faster cost-effective 
regular fixed-point microprocessors. The new approach to the 
design is based on a technique of numerical errors propagation 
suppression applied to a direct computation method and 
realised by means of controlled massive parallel error- 
correcting feedback loops. 

With regards to stochastic quadratic optimisation, this work 
can be considered as an approach to filling a gap between 
slow-convergent but stable gradient-based searching 
algorithms featured with parallel error-correcting feedback 
loops and fast-convergent but numerically instable direct 
methods, which are known to be without effective 
mechanisms of feedback error-correction. 

The rest of the paper is organized as follows. Section II 
presents the conjugate direction decomposition (CDD) of the 
optimum filter weights. In Section III, we derive a new 
stochastic fast-convergent numerically stable STAP algorithm 
that computes the metric-orthogonal conjugate-direction basis 
using a new version of the modified Gram-Schmidt 
orthogonalization (MGSO) method. Results of numerical 
simulations for several sizes of fixed-point binary 
data/operation representation are demonstrated in Section IV. 
Brief conclusions are summarized in Section V. 

II. Background 

In order to maximize the SINR at the space-time filter 
output for a given amount of input samples, we have to solve 
the generalized least-squares (GLS) problem [8] 

rninll v 5 -R^w„ || 2 , n=\,2,... (i) 

where w n is the A^-vector of filter weights, v s is the reference 
A^-vector and R Xn is an estimator of the interference-plus-noise 
covariance matrix obtained using n input sampled /^-vectors x,. 

Using the singular value decomposition (SVD), the GLS 
solution can be represented as [8] 

w ; =r;„v s =u a:u> s =Z^v s |> " ( 2 ) 

where Ry„ + is the Moore -Penrose pseudoinverse K*K matrix, 
U„ is the K*r n matrix of singular vectors and A„ is the real 



diagonal matrix of the corresponding singular values arranged 
in descending order. The number of non-zero singular values, 
r„ - the effective rank of matrix R x „ - is always bounded, i.e. 



r <mm(n,K) 



(3) 



According to the solution (2), the space-time adaptive filter 
is a two-stage A^-input signal processor that contains a 
decorrelation filter associated with matrix U„ followed by a 
.K-input optimum linear combiner. 

The STAP algorithms developed using the EVD/SVD 
approach are known to be unsuitable for implementation on 
real-time fixed-point processors due to the requirements of 
arithmetic precision and numerical complexity [9-11]. 

In this paper, we suggest a new factorisation of the 
optimum GLS solution using an alternative vector basis that 
allows us to derive a new family of fast convergent adaptive 
algorithms with improved numerical stability, better fault 
tolerance, and lowered computational complexity. 

Linear algebra states [12] that the optimum weights in (2) 
can also be represented by making use of a basis of conjugate 
directions as 



w + =S A + S H v t 



-l.fi. 



(4) 



where S„ is a Kxr„ column matrix containing r n conjugate 
directions which are orthogonal only with respect to the 
matrix R Xn , i.e. 



and 



s„, R a-„s"=0, i,je[l,rj 



A„ =ting{SJ! =diag{s"R A ,,s„ i K' 



(5) 
(6) 



contains the R^-norms of conjugate directions arranged in 
diminishing order. 

From (5) it is evident that the matrix S„ decorrelates the 
filter input signals as well as the matrix U„ thus making the 
elements of the output vector y n = S"x n mutually orthogonal, 

i.e. R K =£{y n yf}=A„. 

Although the conjugate directions are orthogonal only in 
the metric described by the input covariance matrix, the CDD 
offers one clear advantage over the EVD/SVD. Given a KxK 
nonnegative-definite matrix R x „ and an arbitrary set of K 
linearly-independent A'-vectors, the conjugate basis can be 
computed in K recursions utilizing simple stable inexpensive 
methods of vector orthogonalization. 

III. Adaptive Algorithm Based on the Modified Gram- 
Schmidt Orthogonalization 

In the first part of this section, we introduce a new version 
of the modified Gram-Schmidt orthogonalization method 
(MGSO) to compute a metric -orthogonal vector basis with 
suppressed levels of errors accumulation and propagation. A 
novel fast stable STAP algorithm developed basing on the 
principle massive staggered feedback orthogonalization is 
presented in the second part. 



A. Numerically Stable Metric-Orthogonalization Method 

Consider a KxM column matrix S containing an arbitrary 
set of M linear-independent vectors (M < K) in C K and R is a 
KxK full-rank matrix. Assume also that the vectors s* are 
arranged in order of R-norm decrease. 

In order to construct such a set of R-orthogonal vectors, we 
use the following recursive orthogonalization form that differs 
from the MGSO only within the metric of orthogonality 



S(i) = S(i-1)H(i), *' = 1,...,M-1 



(7) 



where S(i) is the KxM matrix containing i+l conjugate 
directions and H(7) is the MxM basis transformation matrix, 
which differs from the MxM identity matrix only in the z-th 
row whose non-zero elements are calculated as 



sfRs t 
Vrs, 



k = i + l,...,M 



(8) 



Similarly to the stabilized MGSO [13] the remaining 
vectors s, +1 , ... , s M must be rearranged in order of R-norm 
decrease to suppress the volumes of transferred numerical 
errors implying that ||/z ;/ t(0l| 2 < 1 for all i and k. 

Thus the resulting set of R-orthogonal vectors is 



S = S f[H(0 = S H 



(9) 



where H is the MxM upper triangle matrix whose off-diagonal 
non-zero elements are limited in absolute value to unity. 

From (9) it follows that if the initial vector system S is set 
to be the identity matrix columns in the case of M= A'then the 
final conjugate basis becomes S = H having only K(K+l)/2 
non-zero components. 

B. CDD Based Adaptive Algorithm 

In most practical cases, the input covariance matrix R x „ is 
estimated using the following well-known recurrent form 



R, 



: (1-X)R X( „_ 1) +Xx„x:, R xo =0, 



(10) 



According to the recurrent update in (10), the conjugate- 
direction basis for the decomposition in (3) is built using the 
stable version of metric-orthogonalization given above as 



S, =(l-/„)S„_ 1+ XS n H , S =I 4 



(11) 



The structure of elementary R x „-orthogonalization matrices 
H„ in (11) depends on the sample number n and estimated 
effective rank r„. For example, matrix Hi may have only one 
row containing K-l non-zero elements in addition to its unit 
diagonal while matrix H K may contain r„ 'non-zero' rows. 

The elements of an i-th non-zero row in matrix H„ are 



° ' ni 

h nlk =0, if a\ < a\ k , k = l,..., r 



(12) 



where 



= 0- 



S X 



, k = 1, .... 



(13) 



In (13), «/ =2,3, ... is the amount of sequential input 
samples that may be required to smoothen the estimates of the 
output second-order statistics at each recursion, y ni is the ;-th 
output of the matrix S„ defined decorrelation filter, c„a is the 
averaged covariance between the i-th and k-th filter output and 
<j ni is the i-th filter output averaged power. 

Finally, the adaptive MSINR weights in (4) become 
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(14) 



Equations (11) to (14) highlight how the use of the CDD 
has allowed us to replace direct metric-orthogonalization of 
input vectors x„ with regular orthogonalization applied to the 
decorrelation filter outputs y„. In thus way there is effective 
error control on the performance of adaptive space-time 
processing. 

A computation schematic that summarises the design of the 
new fast stable STAP algorithm is presented in Table I. 

TABLE I 

CDD Based Fast Stable STAP Algorithm 



♦ Initialisation 

- transformation matrix: H = l K 



- matrix of conjugate directions: S : 
♦ Adaptation 
for n = 1,2, ... do 

- current amount of the CD's: N s = n; if n > K then N s = K 

- initialized covariation matrix at the decorrelator's outputs: 
Rro = ®k 

- smoothing the decorrelator's output covariance matrix: 
for i= 1, do 

R yn = (1 - X)R,,„_ 1) +XSl 1 x„ H)t / (< _ 1H ,S„ 4 

end i 

- finding the current index to the next maximum output power: 
isn = argmax {r m } , k = 1 , . . . , K, k + i Sn . h i Sn _ 2 , . . . 

- updating matrix H„: 
for m = i su ... , i s „Ao 

for k = 1 , . . . , K do 

if r Ymm > rjkt then h nmk = - r Ynmk lrY nmm else h nmk = 0; 
end k 
end m 

- updating the CD basis: 

S„ =(1-X)S„_,+XS„_ 1 H, 

- updating the weight vector: 



end n 




Fig. 1. Processing structure of space-time adaptive filtering based on 
feedback orthogonalization 



A block diagram of the STAP filter implementing the new 
CDD based algorithm is presented in Fig. 1 . 

From equations (11) to (14) and Fig.l it becomes evident 
that the main feature of the newly designed STAP algorithm is 
in massive feedback orthogonalization. This technique 
effectively controls the depth and stability of adaptive 
decorrelation, and places stringent limits on the decorrelation 
weights thus suppressing accumulation and propagation of 
processing errors and faults to prevent the algorithm against 
explosive divergence. 

IV. Simulations 

Simulations were performed in order to characterize the 
convergence performance of the new MSINR adaptive 
beamforming algorithm in a realistic situation for standard 
formats of fixed-point binary processing arithmetic. 

In the simulations, we examine a uniform half-wavelength 
spaced linear antenna array with K = 16 dipole elements. We 
assume the noise to be spatially isotropic zero-mean white 
Gaussian with unit variance. 

The desired signal direction of arrival (DOA) is (f> s = 0°, its 
power is set to lOdB over the noise floor. All M = 4 
interfering signals are modelled by narrowband zero-mean 
Gaussian random processes, their powers and DOA's are 
presented in Table II. 

TABLE III 



m 


1 


2 


3 


4 


Pm (dB) 


40 


40 


30 


50 




-12 


10 


18 


23 



For all simulation scenarios, the number of iterations was 
set to n F = 100 while the number of integrations required to 
estimate the output statistics in (13) was set to «/ = 2. In order 
to demonstrate the new algorithm's stability of convergence, 
no ensemble averages were supposed. 




The learning curves and antenna beam patterns obtained for 
the 12-bit and 8-bit fixed-point binary data/operation format 
are presented in Fig.2 and Fig.3, respectively. For comparison 
purposes, we also show the SINR learning curves and patterns 
obtained for 16-bit fraction floating-point format (blue color). 
The pre-adaptive beam pattern is marked with black color in 
the b) plots. The vertical red lines in the beam-pattern plots 
indicate the DOA's of the interferences. 

V. Conclusion 

The CDD based approach applied to STAP algorithms 
designs has allowed us to synthesize a new fast-convergent 
adaptive filter with the highest immunity to numerical errors 
and faults that facilitates its real-time implementation on 
substantially faster and cheaper regular fixed-point processors. 

Perfected numerical attributes of the new algorithm are 
provided by its inherent system of feedback orthogonalization 
which controls the performance of space-time decorrelation 
and suppresses propagation of processing errors and faults. 

Due to capability of effective fixed-point implementation, 
this algorithm is expected to outperform most well-known 
STAP algorithms in terms of actual time of convergence. 

To some extent, the newly designed CDD-algorithm can be 
considered as a seamless combination of a fast-convergent 
direct numerical method used to obtain a rough approximation 
to the optimum solution and a searching feedback algorithm 
targeted to improve the achieved approximation. 
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